Temporal response of nonequilibrium correlated electrons 
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Abstract 



In this work we examine the time-resolved, instantaneous current response for the spinless Falicov-Kimball model at 
-half-filling, on both sides of the Mott-Hubbard metal-insulator transition, driven by a strong electric field pump pulse. 
O |The results are obtained using an exact, nonequilibrium, many-body impurity solution specifically designed to treat the 
out-of-equilibrium evolution of electrons in time-dependent fields. We provide a brief introduction to the method and its 
'computational details. We find that the current develops Bloch oscillations, similar to the case of DC driving fields, with 
an additional amplitude modulation, characterized by beats and induced by correlation effects. Correlations primarily 
inanifest themselves through an overall reduction in magnitude and shift in the onset time of the current response with 



1 increasing interaction strength. 
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We examine the time-domain current response for the 
spinless Falicov-Kimball (FK) model at half-filling. This 
■model provides a relatively simple test-bed for understand- 
ing the physics of correlated electron systems. In particu- 
lar, the model possesses a Mott-Hubbard metal-insulator 
transition (MIT) at half-filling, as well as non-Fermi liq- 
uid behavior. [Ij The appeal of this model for the present 
study rests on the exact, nonequilibrium, many-body im- 
purity solution that has been developed for the situation 
where the electrons are subject to time-dependent driving 
•fields. 0^3 Here, the field, or pump pulse, that drives the 
System out of equilibrium is modeled as a large amplitude 
electric field with a characteristic Gaussian profile. In prin- 
ciple the field may take any form, and we also choose to 
study fields modified through the inclusion of a harmonic 
'modulation. This additional modulation leads to notice- 
able changes in the time-resolved current response of the 
"system. 

Nonequilibrium dynamical mean-field theory (DMFT) 
is employed to obtain the real-time dynamics for the FK 
model on the hypercubic lattice in infinite dimensions {d = 
oo). M The Hamiltonian in equilibrium takes the form 



— y 



(c|cj + h.c.) - M ^ c|ci + Uy^Wi c\ci 



The FK model describes the hopping of conduction elec- 
trons that experience an on-site interaction U with an- 
other species of localized electrons where Uc — ^/2t* is 



the critical interaction strength for the MIT at half-filling. 
Throughout this work, the energy unit is taken to be t* . 
At half-filling, the average localized electron occupation 
{wi) — 0.5 and the chemical potential /i = U/2. The 
nonequilibrium DMFT formalism proceeds in essentially 
the same manner as the iterative approach applied in equi- 
librium where all quantities now have two time indices. 

The driving term is modeled by a spatially uniform 
electric field along the (1,1,1,...) hypercubic body diago- 
nal. This choice simplifies evaluation of the noninteracting 
Green's function (GF). [i,[3| The field takes the form 



E(i) = Eo cos(w(< - to) + (j>) exp' 



The influence of this driving field enters through the as- 
sociated vector potential A(t) in the Hamiltonian gauge 
using a Peierls' substitution [6|. The time- resolved current 
response of the system is the velocity weighted average of 
the equal-time lesser GF 0, Q 
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where all velocity components Vi are equal for the cho- 
sen driving field direction. In the nonequilibrium DMFT 
formalism the sum over momentum can be converted to 
a two dimensional integral over two "band energies" dis- 
tributed with a Gaussian joint density of states on the 
hypercubic lattice. The velocity components also can be 
expressed in terms of these band energies. Q The con- 
tinuous matrix operators of the nonequilibrium DMFT 
formalism are approximated by discretizing the Keldysh 
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contour. Quadratic extrapolation of results to zero "step 
size" on the contour ensures that sum rules for the spectral 
moments are satisfied within a few percent for strong cor- 
relations and high fields. 0] This quadratic extrapolation 
based on Lagrange interpolation also produces accurate 
results for current. 

The Gaussian distribution in both band energies makes 
Gaussian quadrature an ideal method for dealing with the 
two dimensional Hilbert transform that yields the local 
GF. In this study the results for the local GF and local 
self-energy were obtained by averaging the Hilbert trans- 
form over quadratures with TV = 54 and = 55 points 
in each dimension for a total of 5,941 quadrature points. 
This represents the most computationally intensive por- 
tion of the algorithm, but also provides an opportunity 
for parallelization as the evaluation of the integrand of the 
transform is independent for every quadrature point. 

In practice, the integrand for each quadrature point, 
obtained from a matrix inversion and two matrix mul- 
tiplies, is a dense, complex matrix whose size depends 
on the discretization of the Keldysh contour. For the 
cases considered here, reasonable accuracy for the spec- 
tral moments following quadratic extrapolation could be 
achieved using matrices ranging in size up to 4100 x 4100 
elements. [Sj] This provides a rather severe "bottleneck" to 
efficient scaling, using the master-slave model, if each slave 
were to communicate its results directly with the master. 
Instead, a recursive binary gather 8] has been employed 
which involves a progressive binary division of the slaves 
and inter-slave communication to accumulate the results 
of the quadrature integration for the Hilbert transform. 
Additional details about the numerical issues associated 
with this method may be found in Refs. [8l-[lo|. Note that 
in this form the parallel code scales reasonably well on a 
number of different platforms with the limitation 

coming from the serial portion of the code where the im- 
purity solver determines the local, double-time self-energy. 
The results presented from this study have been obtained 
from calculations performed on typically ^ 1000 — 1500 
processors with ^ 30 — 40 nonequilibrium DMFT itera- 
tions to achieve convergence of the self-energy. 

Figure [T] shows the time-resolved current response for 
the FK model at half-filling, and for various interaction 
strengths, subjected to a pump pulse with parameters 
Eo = 20, to = 0, cr = 5, w = 0, and = 0. Bloch oscil- 
lations are visible clearly in each panel, especially within 
the insets showing the current response near the center of 
the pump pulse at t = 0. The frequency of these Bloch 
oscillations, proportional to the electric field strength E 
for DC driving fields, 0, 0] varies as a function of time 
due to changes in the field strength associated with the 
Gaussian profile of the pump pulse. An additional ampli- 
tude modulation, characterized by beats in the current re- 
sponse, also appears with a frequency proportional to the 
interaction strength U . The Bloch oscillations and am- 
plitude modulation appear more regular for interactions 
strengths U > Uc than for the case of a DC driving field. 
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Figure 1: Time-resolved current response of the FK model with in- 
teraction strengths (a) U = 0.5, (b) U = 1.0, (c) U = 1.5, and (d) 
U = 2.0 for a system subject to an electric field pulse whose profile, 
plotted in arbitrary units, is shown in panel (a) {Eo = 20, to = 0, 
a- = 5,a; = 0, <j) = 0). The inset in each panel highlights the current 
response near the maximum in the field profile at a time delay t = 
(between time delays t = —2.5 and 2.5). (Color online.) 
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Figure 2: Time-resolved current response of the FK model with in- 
teraction strengths (a) U = 0.5 and (b) U = 2.0. The profile of the 
electric field pulse, plotted in arbitrary units, is shown in panel (a) 
{Eo = 20, to = 0, a = 5, uj = 0.5, </> = 0). The inset in the two 
panels highlights the current response between time delays t = —2.5 
and 2.5. (Color online.) 



This behavior hkely resuhs from the large amphtude of 
the Gaussian pump pulse that reduces the effective damp- 
ing from electron correlations, especially near the center 
of the pump pulse. However, damping effects are manifest 
in the general reduction in the magnitude of the current 
response and a shift in the onset of the current response 
to higher time delays with increasing correlation strength, 
progressing from Fig. [TJa) to Fig. [Hd). 

The results in Fig. [2] show the response for a pump 
pulse with parameters Eq — 20, to — 0, a — 5, u — 0.5, 
and = 0. The pump pulses in this study are assumed 
to be phase-locked (the phase does not vary from pulse 
to pulse or this could be viewed as the "single-shot" re- 
sponse). While one clearly sees Bloch oscillations, they 
appear far more irregular than those presented in Fig. [1] 
This is presumably due to the additional harmonic mod- 
ulation of the pump pulse, especially for U > Uc (see 
Fig. IHb)). For weak correlations, the onset of the current 
occurs in a region where field is strong enough, although 
still quite small compared to the maximum at t = 0. The 
additional amplitude modulation appears for U < Uc in 
Fig. IHa), but this behavior is apparent only within the 
relatively narrow central peak of the pump pulse. The 
damping leads to a more irregular evolution of the Bloch 



oscillations with increasing correlation strength compared 
to the results in Fig. [T] however, there remains a general 
reduction in the magnitude and onset of the current re- 
sponse with increased correlations, even more pronounced 
than that shown in Fig. [TJ 

The results presented from the current model capture 
the formation of Bloch oscillations in the time-resolved 
current response of a correlated system subject to a par- 
ticularly strong pump pulse with a Gaussian profile. Of 
course, this behavior is quite similar to the results pre- 
sented for DC driving fields, 0, Q but the field strength 
is greatly increased to ease issues associated with con- 
vergence and extrapolation in the present study. The 
time-resolved current is one of the most readily accessi- 
ble response functions that can be used to characterize 
the out-of-equilibrium behavior of a correlated electron 
system subject to time-dependent fields. Using the same 
data generated from these calculations, one also can ac- 
cess the time-resolved photoemission response ll|, |l2| that 
provides a more direct measure of the redistribution of 
spectral weight and rearrangement of electronic states that 
should characterize the nonequilibrium correlated system 
in a pump-probe experiment. 
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